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BOUNDED SOLUTIONS TO THE ALLEN-CAHN EQUATION 
WITH LEVEL SETS OF ANY COMPACT TOPOLOGY 

ALBERTO ENCISO AND DANIEL PERALTA-SALAS 


Abstract. We make use of the flexibility of infinite-index solutions to the 
Allen—Cahn equation to show that, given any compact hypersurface E of R d , 
with d ^ 4, there is a bounded entire solution of the Allen—Cahn equation 
on R“ whose zero level set has a connected component difieomorphic (and 
arbitrarily close) to a rescaling of E. More generally, we prove the existence of 
solutions with a finite number of compact connected components of prescribed 
topology in their zero level sets. 


1. Introduction 

The study of the analogies between the level sets of the solutions to the Allen 
Cahn equation 

A u + u — u 3 = 0 

in and minimal hypersurfaces in was greatly fostered by De Giorgi’s 1978 
conjecture that all the level sets of any entire solution to the Allen-Cahn equation 
that is monotone in one direction have to be hyperplanes for d ^ 8. This is a natural 
counterpart of the Bernstein problem for minimal hypersurfaces, which asserts that 
any minimal graph in must be a hyperplane provided that d ^ 8. Ghoussoub- 
Ghi and Ambrosio-Cabre proved De Giorgi’s conjecture for d = 2,3 [Hi, and 
the work of Savin m showed that it is also true for 4 ^ d ^ 8 under a weak 
additional technical assumption. Del Pino, Kowalczyk and Wei [7] employed the 
Bombieri De Giorgi-Giusti hypersurface to show that the statement of De Giorgi’s 
conjeture does not hold for d ^ 9. 

In dimension 2, it is well known [6] that the monotonicity hypothesis can be 
relaxed to the assumption that the solution u is stable, i.e., that its Morse index 
is 0. Let us recall that the Morse index of u is the maximal dimension of a vector 
space V C C£°(R d ) such that 



v 2 + 3u 2 v 2 ) dx < 0 


for all nonzero v £ V. Remarkably, it has been shown recently m that in dimen¬ 
sion 8 (actually, in any even dimension d ^ 8) there are bounded stable solutions 
to the Allen-Cahn equation whose level sets are not hyperplanes, but rather they 
are asymptotic to a minimal cone. For the role of minimal cones in the Allen-Cahn 
equation, see also [4] and references therein. In dimensions d St 7, the level sets 
of stable solutions to the Allen-Cahn equation are conjectured to be all hyper¬ 
planes m ■ 
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The analysis and possible classification of bounded entire solutions to the Allen 
Cahn equation is an important open problem where the Morse index of the solutions 
plays a key role. Unlike the stable case [6], the structure of solutions with finite 
Morse index can be very complex; in fact, in dimension 3 a result of del Pino, 
Kowalczyk and Wei [9] ensures that, under mild technical assumptions, given any 
embedded complete minimal surface in R 3 with finite total curvature, there is a 
bounded entire solution to the Allen-Cahn equation with a level set that is close to 
a large rescaling of this minimal surface, and that the Morse index of this solution 
coincides with the genus of the surface. Also in this direction, the existence of solu¬ 
tions to the Allen-Cahn equation with a level set close to a nondegenerate minimal 
hypersurface was proved by Pacard and Ritore [14] provided that the ambient space 
is a compact Riemannian manifold (instead of R d ). Furthermore, Agudelo, del Pino 
and Wei [2] have recently constructed bounded entire axisymmetric solutions on R 3 
of arbitrarily large index that have multiple catenoidal ends. 

Generally speaking, it is expected M that the condition that the Morse index 
of the solution be finite should play a similar role as the finite total curvature as¬ 
sumption in the study of minimal hypersurfaces in Euclidean spaces. In particular, 
it is well known that there are many infinite-index solutions to the Allen-Cahn 
equation 0®, and this abundance of solutions should translate into a wealth of 
possible level sets. 

Our objective in this paper is to explore the flexibility of bounded entire solutions 
to the Allen-Cahn equation of infinite index by showing that there are bounded 
solutions to the Allen-Cahn equation on R d with level sets of any compact topology. 
Specifically, given a compact hypersurface E without boundary of R d , we will show 
that there is a rescaling of E that is arbitrarily close to a connected component of the 
nodal set of a bounded entire solution of the Allen-Cahn equation. Furthermore, 
this level set is structurally stable in the sense that any function on R d which is 
sufficiently close to u in the C 1 norm in a neighborhood of this set will also have 
a zero level set of the same topology. In view of the existing literature, we are 
particularly interested in the case of high dimension d. 

To present a precise statement, let us agree to say that an e-rescaling is a diffeo- 
morphism of R d that can be written as $ = $1 o $ 2 , where $2 is a rescaling and 
||$i — id|| c i( R d) < e (here we could have taken any other fixed C k norm, though). 
By a hypersurface we will refer to a smoothly embedded codimension 1 submanifold 
of R d , so self-intersections will not be allowed. Furthermore, in what follows we 
will use the notation (a;) := (1 + l^l 2 ) 1 / 2 for the Japanese bracket. 

Theorem 1.1. Let E be any compact orientable hypersurface without boundary 
of R d , with d ^ 4, and take any e > 0. Then there is an entire solution u of 
the Allen-Cahn equation in R d such that its zero level set u -1 (0) has a connected 
component given by 4>(E), where $ is an e-rescaling. Furthermore, this set is struc¬ 
turally stable and u falls off at infinity as |u(x)| < C{x)~*~. 

It is worth mentioning that the result that we will actually prove (Theorem 14.II) 
is in fact stronger, in the sense that given any finite number of hypersurfaces 
Ei,..., Ejv that are not linked (see Definition 12.111 we will show that there is a 
diffeomorphism $ such that 4>(Ei) U • • • U $(Ejv) is a union of connected compo¬ 
nents of the nodal set of a bounded entire solution to the Allen-Cahn equation. The 
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diffeomorphism $ is not an e-rescaling, although it does act on each hypersurface Ej 
as an e-rescaling composed with a rigid motion. 

The idea of the proof of the theorem is that, when u is small in a suitable sense, 
solutions to the Allen-Cahn equation behave as solutions to the Helmholtz equation 

Aw + w = 0. 

Hence a key step of the proof is to establish an analog of Theorem 1 1.1 1 for solutions 
to the Helmholtz equations with the sharp fall-off rate at infinity, which is as ( x ) t - 
fTheorem 12.21) . For this we combine a construction using the first eigenfunction of 
the domain bounded by E with a Runge-type theorem with decay conditions at 
infinity that generalizes the results that we proved in nans] for Beltrami fields 
on R 3 . Using suitable weighted estimates for a convolution operator associated 
with the Helmholtz equation (Theorem 13.11) , we then promote these solutions of 
the Helmholtz equation to solutions of the Allen Cahn equation and show that the 
latter still possess a nodal set of the desired topology. From the method of proof 
it stems that the statement of Theorem o remains valid for much more general 
nonlinearities. (More precisely, one can replace u 3 by a smooth enough function 
F(u) that behaves as u 1+a as u —> 0 for some a > 0. The statement then remains 
valid provided the dimension is larger than some explicit constant do (a).) 


2. Bounded solutions to the Helmholtz equation 


In this section we will prove an analog of Theorem 11.11 for solutions to the 
Helmholtz equation on R d . We shall begin by introducing some notation. 

Let us consider the function 

(2.1) G{x) := /3 |x| 1- ^ Yd_ x {\x \), 

where Yj__, denotes the Bessel function of the second kind and we have set 

o 1 


P '■= 


i—ii 


r(f-i) 


with |S d-1 | the area of the unit (d— l)-sphere and F the Gamma function. A simple 
computation in spherical coordinates shows that AG + G = 0 everywhere but at 
the origin and the asymptotics for Bessel functions shows that 


G{x) = 



o(M 3 - d ) 


as x —> 0. It then follows that G is a fundamental solution for the Helmholtz 
equation, so if v is, say, a Schwartz function on R d the convolution G * v satisfies 


( 2 . 2 ) 


A(G *v) + G*v = v. 


As we discussed in the Introduction, we will prove a result that is considerably 
more general than Theorem 11.11 as it applies to an arbitrary number of hypersur¬ 
faces. There is, however, a topological condition that we must impose on these 
hypersurfaces, which is described in the following 

Definition 2.1. Let Ei,..., Ejv be compact orientable hypersurfaces without bound¬ 
ary of R d . We will say that they are not linked if there are N pairwise disjoint 
contractible sets Si,..., Sn such that each hypersurface Ej is contained in Sj. 
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We are now ready to state and prove the main result of this section. Notice 
that the proof of the theorem provides a satisfactory description of the structure 
of the diffeomorphism \I>, as noted in Remark 12.31 below. The proof makes use of 
some techniques we introduced in m to study the level sets of harmonic functions 
and in m to construct Beltrami fields with prescribed vortex tubes. Throughout, 
diffeomorphisms are assumed to be of class C°° and connected with the identity, and 
Br denotes the ball centered at the origin of radius R. Observe that, of course, for 
TV = 1 the condition that the hypersurface be not linked is empty, as it is satisfied 
trivially. 

Theorem 2.2. Let Ei,..., Eat be compact orientable hypersurfaces without bound¬ 
ary ofW 1 that are not linked, with d 3. Then there is a function w satisfying the 
Helmholtz equation 

Aw + w = 0 

in R d and a diffeomorphism US' ofM. d such that 'F(Ei),..., ’F(EAr) are structurally 
stable connected components of the zero set u'” 1 (0). Furthermore, w falls off at 
infinity as \d a w(x)\ < C a (x)^~ for any multiindex a. 

Proof. An easy application of Whitney’s approximation theorem ensures that, by 
perturbing the hypersurfaces a little if necessary, we can assume that £ ? is a real 
analytic hypersurface of R d . The fact that the hypersurfaces are not linked allow 
us now to rescale and translate them so that the (unique) precompact domains Llj 
that are bounded by each rescaled and translated real-analytic hypersurface, which 
we will call £' := dflj, are pairwise disjoint and their first Dirichlet eigenvalue 
Ai(flj) is 1. The first eigenvalue is always simple, so there is a unique eigenfunction 
. modulo a multiplicative constant, that satisfies the eigenvalue equation 

Aipj + if j = 0 in flj , ifj | s' = 0 . 

We can choose ifj so that it is positive in flj. 

Hopf’s boundary point lemma shows that the gradient of ifj does not vanish 
on E j\ 

(2.3) min |V0 j(t)| > 0. 

x£Y,j 

Furthermore, as the hypersurface Ej is analytic, it is standard that ipj is analytic 
in an open neighborhood Qj of the closure of f lj. 

Our goal is to construct a solution w of the Helmholtz equation in that 
approximates each function ifj in the set f lj. To this end, let us take a smooth 
function % : —>• M that is equal to 1 in a narrow neighborhood of the closure Q 

and is identically zero outside H, with 

N N 

(1 := (J n,j , n := (J o ; . 

i=i 3 =1 

We can now define a smooth function w i on R. d by setting 

N 

Wi 

3 =1 

Here we are assuming that iv\ := 0 outside fh 
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Since w\ is compactly supported, we can employ the fundamental solution EH 
to write 

(2.4) wi{x) = [ G(x-y) f{y)dy 

Jw 1 


with / := Awi + w±. The support of the function / is obviously contained in 
the open set fi\fi. Therefore, an easy continuity argument ensures that one can 
approximate the integral (12.41) uniformly in the compact set fi by a finite Riemann 
sum of the form 

M 

(2.5) w 2 (x) '■= Y, c n G{x - x n ). 

n= 1 

Specifically, it is standard that for any 6 > 0 there is a large integer M, real 
numbers c n and points x n G fi\fi such that the finite sum (12.51) satisfies 

(2.6) \\wi - w 2 ||c°(n) < <5- 


Let us now take a large ball Br containing the closure of the set fi. We shall 
next show that there is a finite number of points {x' n }^ =1 in R d \BR and constants 
c' n such that the finite linear combination 

M' 

(2.7) W 3 {x) := c n G ( x - x 'n) 

71=1 

approximates the function w 2 uniformly in fh 

(2.8) 11^2 — W 3 ||c°(f2) < S. 

Here 5 is the same arbitrarily small constant as above. 

Consider the space V of all finite linear combinations of the form EH where x' n 
can be any point in M d \Bn and the constants c' n take arbitrary values. Restricting 
these functions to the set H, V can be regarded as a subspace of the Banach space 
C°(fi) of continuous functions on Q. 

By the Riesz-Markov theorem, the dual of C°(fl) is the space A4(fl) of the finite 
signed Borel measures on whose support is contained in the set fi. Let us take 
any measure /j G A4(H) such that f Rd f dy, = 0 for all / G V. We now define a 
function F G Ll oc (R d ) as 

F(x) := f G(x — x') dy(x '), 

J s. d 

so that F satisfies the equation 

AF + F = y . 

Notice that F is identically zero on R d \Bn by the definition of the measure y and 
that F satisfies the elliptic equation 

AF + F = 0 

in R d \fl, so F is analytic in this set. Hence, since R d \H is connected and contains 
the set R d \Bu, by analyticity the function F must vanish on the complement of 
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Cl. It then follows that the measure /r also annihilates the function G(- 
y ^ Q because 

0 = F(y) = [ G(y - x') dy(x'). 

J R d 


Therefore 



W 2 dy, = 0 , 


y) with 


which implies that W 2 can be uniformly approximated on fI by elements of the 
subspace V, due to the Hahn-Banach theorem. Accordingly, there is a finite set of 
points {x' n }^L. 1 in R. d \Bn and reals c' n such that the function w 3 defined by (EH 
satisfies the estimate (12.81) . 


To complete the proof of the theorem, notice that the function u> 3 satisfies the 
equation 


(2.9) 


Aw3 + W3 = 0 


in the ball Br, whose interior contains Cl. Let us take hyperspherical coordinates 
r := |ie| and ui := x/\x\ £ in Br. Expanding the function W 3 (with respect 
to the angular variables) in a series of spherical harmonics and using Eq. (12.91) , we 
immediately obtain that W 3 can be written in the ball as a Fourier-Bessel series of 
the form 

OO 

w 3 = ^2^2 C lmjl( r ) Ylm{u) , 

i=0 m£l[ 

where ji denotes a d-dimensional hyperspherical Bessel function, Yi m are spherical 
harmonics on and Ii is a finite set that depends on l and whose explicit 

expression will not be needed here. 

Since the above series converges in L 2 (Br), for any 5 > 0 there is an integer l 0 
such that the finite sum 

lo 

w := ^2 ^2 Y i™{u) 

1=0 

approximates the function W 3 in an L 2 sense: 

( 2 . 10 ) \\w - w 3 || l 2 ( Br ) < S. 

By the properties of Bessel functions, w is smooth in K. d , falls off as 

\d a w{x)\ < C a (x) L Y i 

at infinity for any multiindex a and satisfies the equation 

(2.11) Aw + w = 0 
in the whole space. 

Given any R' < R large enough for the set Cl to be contained in the ball Br', 
standard elliptic estimates allow us to pass from the L 2 bound (12.101) to a uniform 
estimate 

\\w - w 3 || co ( b h ,) < C6. 

From this inequality and the bounds (12.61) and (12.81) we infer 

(2.12) ||w — wrllco(n) < C6 . 
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Moreover, since wi also satisfies the Helmholtz equation in a neighborhood of 
the compact set fi, standard elliptic estimates again imply that the uniform es¬ 
timate (12.121) can be promoted to the C 1 bound 

(2.13) ||u> - Willcqn) < CS. 

Finally, since £i U • • ■ U Sat is a union of components of the the nodal set of 
wi and the gradient of wi does not vanish on these hypersurfaces by (12.311 . the 
estimate (12.131) and a direct application of Thom’s isotopy theorem [T] Theorem 
20.2] imply that there is a diffeomorphism if of M. d such that 

(2.14) T(SiU---US A r) 

is a union of components of the zero set w _1 (0). Moreover, the diffeomorphism T 
is C 1 -close to the identity. The structural stability of the set (12.141) for the function 
w also follows from Thom’s isotopy theorem and the lower bound 

min | Vw(x)\ > 0 , 

xe^EiU—UEiv) 

as a consequence of the G 1 estimate (12.131) and the fact that the function w\ satisfies 
the gradient condition (12.31) . □ 

Remark 2.3. It follows from the proof that there are rescalings \E'|, translations 
and diffeomorphisms with ||Tj — id||< 7 i( R <q arbitrarily small such that 

$(E i ) = (^o^ 0 $3 )(Sj . ) _ 

In particular, if N = 1 the diffeomorphism T can be assumed to be an e-rescaling. A 
minor modification of the argument would have allowed us to take ||^/j — id||c 7 *(H. d ) 
arbitrarily small, with k any fixed number. 

3. A WEIGHTED ESTIMATE FOR A CONVOLUTION OPERATOR 

In promoting solutions to the Helmholtz equation with sharp decay at infinity 
to solutions to the Allen-Cahn equation, the estimates that we establish in this 
section will play a key role. 

Specifically, we will be interested in the convolution of the fundamental solu¬ 
tion G, introduced in Eq. m, with functions with certain decay rate at infinity. 
To quantify this, for any nonnegative integer k and any positive real v let us denote 
by G^(R d ) the closure of the space of Schwartz functions on with respect to the 
metric 

IMU,„ := max sup \(x) v d a v(x)\. 

Clearly 

\\vw\\k, v + v > < C'NMHk*'/ 

whenever v G C*( R d ) and w £ C^, (R d ), where C is a constant that only depends 
on k. In particular, 

(3.1) IKIU,,<C||u||^ /s . 

The following theorem, which asserts that the convolution with G defines a 

bounded map G^(R d ) —> C l ^_ 1 (R d ) for any v > d, provides the estimates that we 

2 

need: 
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Theorem 3.1. Suppose that d ^ 3. Then for any v £ C k (M. d ) with k ^ 0 and 
v > d, one has 

||G * v\\ k d^i < CIMU^ 

with a constant that depends on d and v but not on v nor k. 


Proof. In view of the well known asymptotics for Bessel functions when d ^ 3, 
there is a positive constant C such that G is bounded by 


It then follows that 


|G(x)| < 


C\x\ 2 ~ d 
C |x|^ 


if |x| < 1, 
if |x| > 1. 


|G*x(x)|< f \G{z)\\v{x — z)\dz 
JwL d 


(3.2) 


<G|M|o,„ 



z) u dz+ f | x | 12 d 
J& d \B 1 



For any fixed x, the first integral is convergent for any value of u, while the second 
converges provided that v > Since v > d > we infer that G * v(x) is well 
defined as a convergent integral for any v € C°(M“) and all x £ R d , and it only 
remains to analyze its behavior for large |x|. 

For concreteness, let us assume that |x| > 2. We shall next show that the 
integrals 


h := f \z\ 2 ~ d (x - z)~ v dz 
Jb x 

h '■= f \z\~^~ {x — z)~ v dz , 

dB \ x \/2 

I 3 := f \z\^ (x-z)-^dz, 

J B 2 \x\ \-S|a: 1/2 

I 4 '■= f \z\~ r {x — z)~ v dz 
JM. d \B 2 \ x \ 

are then bounded as 

(3.3) Ij < G|x|^ , 

where C does not depend on v. In view of the inequality (13.21) and the fact that 


1 -d 
2 


(x 


z) V ^ I 2 + ^3 + I 4 , 


this shows that the convolution with G is a bounded map G°(R d ) —» G°_i(R d ). 

~T~ 

Since for any multiindex a we have 


<9“(G * x) = G * (d a v ), 


this immediately implies that the convolution with G is also a bounded map 
G^(M d ) —> C k d _ 1 (R d ), thereby proving the theorem. 
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So it only remains to prove the estimate (13.311 for 1 ^ j ^ 4. For this we will 
start by using the elementary inequality 


( x - z ) > 4 1 


\\ x \ if kl <¥, 


2 \ z \ if \z\ > 2|x| 


to obtain, for Id > 2, 


h < G\x\~ v f \z\ 2 ~ d dz = C\x\~ u < C\x\~ d < G\x\^ , 

J Bi 

I 2 < C\x\~ u f \z\^ dz = C\x\^ L ~' J < C\x\^ , 

J B i^i/o 


dz = <C|®|V 


h<C 

jR d \B 2 \ x \ 

To obtain these bounds we have used that v > d by assumption. 

To estimate I 3 we choose a Cartesian basis so that x = |x| ei and then use the 
rescaled variable z := z/\x\ to write 


h = 


>b 2M \b m/2 


\z\ 2 (x-z) " dz 


(3.4) 


I I d+1 

= \x\ 2 


\z | 2 dz 


>b 2 \b 1/2 (jyp + l e i - A 2 ) 2 


Denoting by B' the ball centered at ei of radius y, one can check that 


III ^ dz 


- < C 


l ‘ 1/4 p d ~ l dp 


b> (|yp + ki -z\ 2 )2 J 0 (j yp+P 2 ) 2 


M 2 

/■M/4 -d- ij- 

= cw l fcw 

poo —d— 1 j — 

< cm-' / -L * 

Jo (1 

< C\x\ v ~ d 


P 2 ) 


where we have defined p := |rr|p, and we have used that the integral in p is conver¬ 
gent for any v > d. Plugging this in (13.411 . one gets that 


h = \x\^r- v 


\z\-t- dz 


Izl^dz 


B' (+ l e i “ z \ 2 ) 2 Jb 2 \(b 1 /2 ub>) (j^p + |ei - z\ 2 ) 2 

<C\x\^ L - u {\x\ v - d + C) 

<C\x\^ . 

To obtain the first inequality we have used that |ei — z \ 2 ^ yg for all z € B 2 \(B 1 / 2 U 
B'), and the second inequality follows from the assumption v > d. This is the last 
estimate that we needed in (13.31) and thus the theorem follows. □ 

In view of the structure of the nonlinearity of the Allen-Cahn equation, the 
following corollary will be useful: 
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Corollary 3.2. For any v £ C^-i (R d ) with k ^ 0 and d ^ 4, one has the estimate 

2 

IIG * < C'lklH^ - 

Proof. We can apply Theorem 13.11 with v := 3rf ~ 3 because u > d for all dimensions 
d ^ 4, thus implying that 

l|G* (t 3 )!!^ < q|u 3 || fei ^3 < c||u|| 3 d_i, 

where we have used the relation m- □ 

4. Proof of Theorem o 

We are now ready to prove the main result of this paper, which reduces to 
Theorem D when N = 1. 

Theorem 4.1. Let Si,..., Ejv be compact orientable hypersurfaces without bound¬ 
ary of W l that are not linked, with d ^ 4, and let us take any positive integer k. 
Then there is a dijfeomorphism $ of R d such that 4>(Ei),..., <f>(£iv) are connected 
components of the level set u -1 (0) of a smooth solution to the Allen-Cahn equation 
in that is bounded as \d a u(x)\ < C a {x)^~ for any multiindex with |a| < k. 
Furthermore, these level sets are structurally stable and the dijfeomorphism can 
be assumed to have the same structure as in Remark \2.H[ 

Proof. By Theorem 12.21 there is a solution w to the Helmholtz equation 

Aw + w = 0 

on R d such that ^(Sj),..., ^(Ejv) are connected components of its zero set w _1 (0), 
where T is a diffeomorphism of R d . Moreover, ||w|L a- 1 < C and the above hyper¬ 
surfaces are structurally stable in the sense that there exist a large ball Br and a 
positive constant such that, if w' is any function with 

(4-1) \\w- w'\\ C i(b r ) <V, 

then there is a diffeomorphism $ of R d such that 

(4.2) $(S 1 )U---U$(E i v) 

are structurally stable connected components of the level set u/ _1 (0). Furthermore, 
$ is close to T in the norm C 1 (R d ). 

Let us take a small positive constant e that will be fixed later and consider the 
iterative scheme 

u o := 6 w , 

(4.3) u n+ i :=6w + G* (u 3 ), 
where we have set 

S:= nw\\ k ^- 

Our goal is to show that if e is small enough, u n converges in C|_i(R d ) to a 

2 

function u that satisfies the Allen-Cahn equation 

A u + u — u 3 = 0 

and is close to Sw in a suitable norm. 
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A first observation is that, if ||u n || fe < e and e is small enough, by the 
definition of S we automatically have 

IK+illfc,^ s Ikllfc^ + ||G * (u 3 n ) || fci d-L 

^ d||w;|| fc d~i + G11 u n 11? d -1 

<l + c ' 3 

(4.4) < e. 


Here we have used Corollarv l3.2l to estimate G* (m^ ). Notice that the smallness that 
we have to impose on e only depends on the constant that appears in Corollarv l3.2l 
In particular, since the first function u o of the iteration satisfies 

ll u o|| fci ^i = 2 > 

the induction property (14.41) then implies that 


(4.5) IKIU^ < e 

for all n. 

To estimate the difference u n +i—u n , let us start by noticing that for any functions 
v. v' we have 


\ k 3d —3 = max sup ( x ) 2 \d a (v 2 (v — v') + vv'(v — v') + v ,2 (v — v')) | 

’ 2 laKfe^gRd 

II 2 


<G(|M|^ + |K||^)||u-u'||,^ i . 


It then follows from Theorem 13. II the fact that :id , 2 3 > d when d ^ 4, and Eq. fi31) 
that we can write 


||Mn+l — IIG * {v? n M^_i)|| fej d^l 

< G||u^ - M^-lllfc, 2 ^ 2 

(4.6) ^ Ge ||u n U n — lH^ AzA ■ 

If e is small enough for Ge 2 < 4, it is standard that m and m imply that u n 

converges in C k d _ 1 (R d ) to some function u with 
2 

(4.7) < e- 

Since the map v K > G * ( v 3 ) is continuous in C k _ i (K. d ), from (14.31) we infer that u 

~T~ 

satisfies the integral equation 

(4.8) u = 6w + G * (u 3 ). 

As ic is a solution of the Helmholtz equation and G is a fundamental solution 
satisfying m , it then follows that 

Am + u = u 3 , 

so m is a solution of the Allen-Cahn equation, which is smooth by elliptic regularity. 
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One can now use the bound (TOT) . the relation (14.811 and the definition of S to 
write 


w — 


= - \\Sw-u\\ k d-i = - ||G * (u 3 ) 


C, 


In view of the stability estimate im if e is small enough (namely, Ce 2 < 77 ), 
we infer that there is a diffeomorphism $, close to the diffeomorphism 4' in the 
norm C 1 (R d ), such that the hypersurfaces (|4.2D are structurally stable connected 
components of the level set zi _ 1 (0). The theorem then follows. □ 


Acknowledgments 

The authors are supported by the ERC Starting Grants 633152 (A.E.) and 335079 
(D.P.-S.). This work is supported in part by the ICMAT-Severo Ochoa grant SEV- 
2011-0087. 


References 

1. R. Abraham and J. Robbin, Transversal mappings and flows, Benjamin, New York, 1967. 

2. O. Agudelo, M. Del Pino, J. Wei, Solutions with multiple catenoidal ends to the Allen—Cahn 
equation in R 3 , J. Math. Pures Appl. 103 (2015) 142-218. 

3. L. Ambrosio, X. Cabre, Entire solutions of semilinear elliptic equations in R 3 and a conjecture 
of De Giorgi, J. Amer. Math. Soc. 13 (2000) 725—739. 

4. X. Cabre, J. Terra, Saddle-shaped solutions of bistable diffusion equations in all of R 2m , J. 
Eur. Math. Soc. 11 (2009) 819-843. 

5. E.N. Dancer, New solutions of equations on R n , Ann. Sc. Norm. Super. Pisa 30 (2001) 535-563. 

6. E.N. Dancer, Stable and finite Morse index solutions on R n or on bounded domains with small 
diffusion, Trans. Amer. Math. Soc. 357 (2004) 1225-1243. 

7. M. del Pino, M. Kowalczyk, J. Wei, On De Giorgi conjecture in dimension n ^ 9, Ann. of 
Math. 174 (2011) 1485-1569. 

8. M. del Pino, M. Kowalczyk, J. Wei, On De Giorgi’s conjecture and beyond, Proc. Natl. Acad. 
Sci. 109 (2012) 6845-6850. 

9. M. del Pino, M. Kowalczyk, J. Wei, Entire solutions of the Allen-Cahn equation and complete 
embedded minimal surfaces of finite total curvature in R 3 , J. Differential Geom. 93 (2013) 
67-131. 

10. A. Enciso, D. Peralta-Salas, Knots and links in steady solutions of the Euler equation, Ann. 
of Math. 175 (2012) 345-367. 

11. A. Enciso, D. Peralta-Salas, Submanifolds that are level sets of solutions to a second-order 
elliptic PDE, Adv. Math. 249 (2013) 204-249. 

12. A. Enciso, D. Peralta-Salas, Existence of knotted vortex tubes in steady Euler flows, Acta 
Math. 214 (2015) 61-134. 

13. N. Ghoussoub, C. Gui, On a conjecture of De Giorgi and some related problems, Math. Ann. 
311 (1998) 481-491. 

14. F. Pacard, M. Ritore, From constant mean curvature hypersurfaces to the gradient theory of 
phase transitions, J. Differential Geom. 64 (2003) 359-423. 

15. F. Pacard, J. Wei, Stable solutions of the Allen-Cahn equation in dimension 8 and minimal 
cones, J. Funct. Anal. 264 (2013) 1131-1167. 

16. O. Savin, Regularity of flat level sets in phase transitions, Ann. of Math. 169 (2009) 41-78. 

Instituto de Ciencias Matematicas, Consejo Superior de Investigaciones Cientificas, 

28049 Madrid, Spain 

E-mail address: aenciso@icmat.es, dperalta@icmat.es 






